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Abstract 

We present an approximation to the matrix element for the process 
e^e~ — > e^e~7, appropriate to the situation where one or both of the 
fermions are scattered over very small angles. The leading terms in 
the situation where all scattering angles are small contains not only 
terms quadratic in the electron mass, but also quartic and even sextic 
terms must be included. Special attention is devoted to the numerical 
stability of the resultant expression. Its relation to several existing 
formulae is discussed. 

For any e^e~ collider, the process of Bhabha scattering is very important. 
Whereas large-angle scattering provides its share of information on QED or 
theories of electroweak interactions, small-angle scattering serves as a tool 
for accurate luminosity determination; in each of these cases, radiative cor- 
rections (and hence bremsstrahlung) cannot be neglected. But the radiative 
process by itself, namely 

e+{pi)e'{qi) ^ e+{p2)e-{q2h{k) , (1) 

is important in its own right, especially when one or either of the electrons 
is scattered over a small, or even negligible, angle. For instance, such events 
form an important background to low-energy photon-photon processes, as 
recently pointed out, for the DA$NE collider, in moreover, interestingly, 
the process (|l|) is, for machines such as LEP, the most important restric- 
tion on the beam lifetimes (see, for instance, M). Since the Bhabha cross 
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section peaks at small fermion scattering angles, and the photon is preferen- 
tially emitted more or less collinearly with a charged particle, such processes 
are dominated by events where all scattering angles are small, and many 
different momentum transfers become small simultaneously. It is therefore 
important to have reliable expressions for the matrix element (we shall always 
understood it to be squared and summed/averaged over all particle spins) 
in these extreme kinematical configurations. An allowed simplification of 
the problem, if we restrict ourselves to small fermion scattering angles, is to 
neglect the s-channel diagrams, and the interference between radiation from 
the positron and that from the electron line. Hence, in what follows we shall 
assume that we only have the two Feynman diagrams in which the photon is 
emitted by the incoming and outgoing positron (radiation from the electron 
is then trivially obtained by e'^ e^). We introduce the following notation: 



p/j. 

s 
t 

Ai,2 

Si 



P-P = 4E^ , 



q-q , 

Pi,2-k , 

iP2 + q2f = {P- kf 



(2) 



The five invariants s, Si, Ai 2, and t, specify the kinematics completely up 
to a trivial azimuthal rotation around the beam axis. Using the algebraic 
program FORM P], we then find for the matrix element, | p: 



AiA2(-t) 



s - Qrn^f + {s + Ani^ - 2A2 + tf 
+ (si - Qm^f + (si + + 2Ai + tf 
- 12m^(t + lOm^) 



5m 



AiAat^ 

2m2 
'Aft^ 



2m^ 



-(Ai-A2)2 



2A2 + t){si - 2m^ - 2Ai + t) 



si - 2m^y + (si - 2m^ + ty + AmH 
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2m^ r n 
-^^[{s-2m^f + {s-2m^ + tf + AmH\ , (3) 

where m is the electron mass. If we request the bremsstrahlung photon to 
have a non-infinitesimal energy, there is only one smallness parameter in our 
problem, namely n = m/E; we shall now discuss what is the order of n of 
the various invariants in the different collinear situations. To this end, note 
that with a relative error of O (yU^) we can write 

\t\= + EE' (1- cos <iquq2)) , (4) 

where E' = E{1 — 2A2/S) is the energy of q2- We can therefore distinguish 
the following collinear situations, in which the cross section displays strong 
peaks: 

. Ci: Ai = 0(/i2), A2,t = 0(/) ; 

. C2: A2 = 0(/i2), Ai,t=O(/i0) ; 

. C3: t = 0(/i2), Ai,2 = 0(/i°) ; 

. C4: Ai,2 = 0(^2),t = 0(/i6) . 

Unless is very close to E, both s and si are O Allowed approximate 
matrix elements are those in which only the leading orders in are kept. 

We can now relate Eq.(H) to the old CALKUL result of |^ by noting that 
in Ci and C2 we may neglect with respect to s, Si, and t, but not to 
Ai 2, so we find 



l-M 



+ sl + z'^ + z\ 



e6 ' ' AiA2(-t) 

2{s\ + z\) IT? 2(5^ + z^ 



z 



Af t2 ^2 ^2 

2A2 - s - t = (pi - q2f + O (/i^) , 
_2Ai-Si-t = (p2-gi)' + 0(yu') , (5) 



which is precisely the CALKUL result for these two diagrams. Obviously, 
this expression is not appropriate to either C3 or C4. We therefore proceed 
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to drop, in Eq.(^), terms that are of the order of m^/AiA2t, which are always 
neghgible, and we find 

(s^ + sl + z^ + zf)/{-/\,/\2t) , 

16m^Z(Ai - A2)/(A?A^t2) 
-16m6(Ai - A^f/iAlAlt^) 
-4m^s,/{Alt) , 
-Am\s/iAlt) , 

-2mVA2 , 

sAi - S1A2 - 2A1A2 . (6) 



Let us tabulate the formal power of /i of each term in the four distinct 
collinear situations, writing ( for the formal power of Z in C4: 





Ai 


A2 


t 


Ao 


Ai 


A2 


A3 




A5 


Ae 


Ar 


^8 


Ci 


2 








-2 


-2 








+2 


-2 


+2 


-2 


+2 


C2 





2 





-2 


-2 








+2 


+2 


-2 


+2 


-2 


C3 








2 


-2 


-2 


-2 





+2 








+2 


+2 


C4 


2 


2 


6 


-10 


2C-18 


-10 


c-14 


-10 


-8 


-8 


-2 


-2 



We note the following points. In order to avoid negative \Ai P, we need to 
have C > 4 in C4: we shall see that this is indeed the case. Furthermore, 
the terms ^5,6,7,8 are always negligible when t is small. This is a reasonable 
approximation once we have assumed that the interference between radiation 
from the e"*" and the e~ can be discarded. Lastly, we have to ensure that 
the form can be evaluated in a numerically stable manner. To this end 
we must control the cancellations inside Z, and in addition prove that the 
matrix element | p is indeed O {fi~^^) in C4. These two last items will be 
solved simultaneously at the end of our discussion. Summarizing the above 
considerations, we propose the following form for the matrix element, valid 



Ao = 

A, = 

A2 = 

A3 = 

A, = 

A, = 

Ae = 

A, = 

As = 

Z = 
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for small t: 



+ AfAft^^^' " - ^^1^2 - 2A1A2) 

(Ai - A2)' . (7) 



Before proceeding, let us compare the expression of Eq.(|^) with existing 
ones. An earlier form of the matrix element, due to Baier et al. and 
Altarelli and coworkers [§, reads 

We see that Eq.(^ differs from Eq.(|^ in that it has a different form for 
Z, and no terms with or — which, however, cannot be neglected, in 
contrast to what is commonly the case with ultrarelativistic matrix elements. 

An alternative approximation to the matrix element was recently pro- 
posed in It includes terms that have an explicit m?, in the following 
combination, for which we use the notation of the present paper: 

^ An? /Si t \2 An? / s t ^ ^ 

Aim? ( 2s{p2 - qif + 3(pi - P2Y 



A1A2 V 

^ s + {P2 - qif - 2{pi - P2f ^ ^ 

Let us first note that this expression is not invariant under time-reversal, 
which in our problem is the operation ^ P2) Qi ^ l2-i ^"^^ ^ —k^. 
What may be even more serious is the peaking behaviour. Expressing G in 
terms of the five independent invariants s, si, t and Ai^2, and retaining only 
those terms that contribute in C3 or C4, we have 

4m2(sAi - siAa)^ 48m2(Ai - Aa)^ 16m2s(Ai - 2m2) 



5 



The first two of tliese terms are O {fJ^'^^) in C4, but tlie last is O since 
in general the combination Ai — 2m^ is O (/z^) rather than O (/i^). We there- 
fore have to conclude that, at least in C4, the result of is not appropriate. 
Note, however, that in the soft-photon limit disucssed in [|], Ai — > so that 
the correct O behaviour is recovered in that case. 

Let us now turn to the numerical stability of the expression (|^). We 
may assume that we can evaluate Ai^2 and t in a stable manner for given 
momenta. It remains to control the cancellation from O (/i^) down to O (fi^) 
inside Z, that occurs in C4. To this end we introduce y = q^/E, u = 1 — y, 
and a vector r'^, by ^2 = "^Qi + ^'^^ Denoting by the subscript L components 
along (fi, and by T the orthogonal components, we find 

r'' = {r^,rT,rL) , 

r° = , 

rL = {-2m^y + t)/{2p) , 

r-r = rn?y'^ + ut , (11) 

with p =|pi |. In C4, then, ty is O (fJ^^), and is O (/i"^), since y = O (/i^). 
Some straightforward algebra now leads to the surprisingly simple form 

Z = s{k-r) . (12) 

The best frame to evaluate this expression is the rest frame of = p2 + k^. 
Denoting the various components in this frame by carets, we have, with 
w = /R^= 0(/i2): 

f° = -{2m'^y + t)/{2w) , 



P = [w"" -m')/{2w) . (13) 
we see that all components of are 0(/x^), and k^ is 0(/i). Since r-r is 
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O (/U^), no further cancellations are to be expected^, and Z is indeed O (/i"^). 



Finally, we prove that, in C4, the matrix element is indeed O In 
this matrix element, the electron line and the virtual photon propagator enter 
via the object 

spins 



^ {q'iq2 + qUi^ + Ig"") ■ (14) 



The last term cancels one of the factors t, and will hence not contribute in 
C4. We now rewrite 2 terms of and r^: 

= (r'^ + g/^)/^ , q^ = (^r^ + uq^yy _ (15) 

Owing to current conservation, the terms with q'^, q'^ evaluate to zero, and 
the relevant part of the electron current is therefore 

- ^ ■ (16) 

This object, when evaluated in the rest frame, is itself of order O 

as we have seen above. In addition, again in the rest frame, the positron 

current is O since there 

Ai,2 = {fi') , 
u{pi,2) = 0(f,'/') . (17) 



This finishes the argument: | ^ is indeed O (fi^^^) in C4, and we can 
evaluate it in a numerically stable manner, without the need for extended 
precision (REAL* 16) as in |^ to control the cancellations. An application of 
the form (|^, to the study of beam-beam-strahlung as a limitation to the 
LEP beam lifetimes, will be presented elsewhere [0. 

ijf -^ag smaller, yet another cancellation might be lurking in the product k-r. 
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